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ACnVE INTERFEROMETRIC SIGNAL ANALYSIS IN SOFTWARE 

Field of the Inventioa 
The invention generally relates to active signal anarlysis systems, and in particular, 
interferometric signal analysis systems. 

Related Applications 

This application is a continuation-in-part of U.S. ^^at App. No. 09/913,132, filed 
Aug. 9, 2001, entitled "Method and System for Signal Detection in Arrayed 
Instrumentation Based on Quantum Resonance", which is a national stage application of 
PCT Application No. USOO/04076 filed on Feb. 17, 2000, Which claims priority from U.S. 
Pat App. No. 09/253,789, now U.S. Pat No. 6,136,541, filed on Feb: 22, 1999, entitled 
"Method and Apparatus for Analyzing Hybridized Biocliip Patterns using Resonance 
Interactions Employing Quantum Expressor Functions, all of which are hereby fiiUy 
incorporated by reference". Also incorporated herein by reference in their entirety are tiie 
following patents based on applications filed contemporaneously with U.S. Pat App. No. 
09/253,789: U.S. Pat No. 6,142,681, entitied: "Mefliod ^d Apparatus for Interpretijig 
Hybridized Bioelectronic DNA Microarray Patterns Using Self-scaling Convergent 
Reverberant Dynamics" and U.S. Pat No. 6,245,511, entitied "Method and Apparatus for 
Exponentially Convergent Therapy Effectiveness MdAitbring Using DNA Microarray 
Based Viral Load Measurements". Also incorporated herein by reference in their entirety 
are the following provisional patent apphcations: U.S. E^at App. No, 60/395,074, filed on 
Jul, 9, 2002, entitled "Active Signal Processing Technique Using Resonance 
Interferometry". U.S. Pat App. No. 60/399,227, filed on Jul. 29, 2002, entitled "High- 
Throughput Quantitation of Signals from Arrayed Platforms Using Resonance 
Interferometry" U.S. Pat. App. No. 60/408,057, filed on Sep. 3, 2002, entitled "Technique 
for Analysis and Quantitation of Arrayed Information Structures". U.S. Pat App. No, 
60/408,844, filed on Sep. 5, 2002, entitied "System ahd Method for Analysis and 
Quantitation of Arrayed Platfonns". U.S. Pat. App. No. 60/440,253, filed on Jan. 15, 2003, 
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entitled "Technique for Detection and Quantitation of Ardiyed Platform Signals". 

Background of the Invention 
1. Overview of Passive and Active Interferometric Signal Analysis 
5 Interferometric analysis has proven useful for performing sigaal analysis, 

particularly for detecting events of interest within signal patterns, such as within a two- 
dimensional image. One example of interferometric analj/'sis involves irradiating a moving 
object, such as an aircraft, with radar pulses, that interact with the object and interfere with 
one anottier and which are then reflected back to a sensing device such as an antemia, 

1 0 Signals received by the antenna are compared against signals originally transmitted in an 
effort to detect a unique spectral signature of the object so as to identify the object This 
form of interferometric analysis is referred to as "active analysis", because signals used to 
irradiate the event of interest interact with the object and with other signals to actively 
modify the resulting interference pattern. In other words, the system being analyzed is 

15 actively modified. In the case of radar, electromagnetic fields surrounding the event of 
interest are actively modified. 

An altemative type of signal analysis is "passive analysis" wherein signals that 
have already been detected are processed in an effort to identify events of interest therein. 
Examples include filtering images in an attempt to enhance features therein. These 

20 techniques are deemed to be passive because no noise of other signals are added into the 
system being analyzed. More specifically, passive interferometric analysis typically 
involves the manipulation of observed intensity patterns or phase information without the 
addition of noise or other signals. . Often, passive analysis systems e?q)loit parametric 
metiiods such as matched filter techniques, maximum a-posteriori (MAP) techniques, 

25 maximum likeUhood estimator (MLE) techniques, singular value decomposition (SVD), 
classical Fourier analysis, parametric distributional clustering (PCA), or median filter 
variants. Non-parametric methods may be employed as well. In addition, quantitation and 
clustering techniques may be employed such as "soft-<jo6iputing" and its variants which 
include neural networks, fuzzy logic systems, and Bayesian inferencing systems. 
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Although a wider range of passive analysis techniques have been exploited, each requires 
that the observed signal have an intensity greater than the intensity of the background 
signal in which it £^pears. In addition, the analysis system must properly estimate the 
backgroimd signal and to properly estimate any variability in the background signal In 
5 some cases, spread spectrum (S-S) amplitude correlation is provided with background 
subtraction, and scaling and normalization is employed to compensate for detector 
nonlinearities. Hence, there are various limitations to passive analysis, particularly tlie 
inability to detect signals having an intensity less than or near to the intensity level of the 
background. 

10 Insofar as the analysis of h5i?ridized biological microarray patterns is concerned, 

heretofore, passive analysis has been the only form of analysis performed. Briefly, a signal 
pattern ou^ut from a biological microarray is normalized to compensate for dye, scaimer, 
or array effects . Then an estimate is made as to the backgroimd signals at feature, 
transcript or chip levels. The background is then subtracted from the foregroimd and a 

15 ratio of foreground minus backgroimd is calculated to quantitate a differential expression 
representative of patterns of interest Often, this calculation is performed using channel 
ratios. In any case, examples of passive analysis of hybridized digitized image patterns 
include PCA, SVD, MLE, hidden Markov Model (HMM), and various statistical 
approaches. As the technique is passive, it is subject to the limitations summarized above, 

20 in particular, only hybridized transcripts with genes above the estimated background can 
be detected. In addition, high caU variability in targets with low expression levels arises 
and high CV in fold change computations occurs. Skewed results can arise when dealing 
with platforms or protocols having a limited dynamic range. Typically, passive 
hybridization analysis techniques require vast amounts of data and exemplars to allow 

25 refinement of statistics driving the call analysis. Heretofore, efforts to improve the 
accuracy and reliability of passive techniques for analyzing hybridized digitized images 
involve wet signal amplification and biological enhancements. Computational 
enhancements of the passive analysis have been fairly minimal. 

Thus, heretofore, computer systems have been largely limited to performing 
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passive analysis whereas active ana])^is is largely performed only by hardware systems. 

2. Hardware-Based Active Signal Analysis 

With hardware-based active interferometric signal analysis, external excitation is 
provided to the system being analyzed while the system is being analyzed to provide an 
5 interaction that generates additional information or provides some change in a qualitative 
property of the system. Typically, the external excitation involves some form of noise 
injection, which results in a nonlinear coiq)ling between the noise and events of interest 
within the system being analyzed, which enhances those events. As noted, one example of 
active signal processing (ASP) is radar wherein the signal being injected into the system 

10 being analyzed is the radar signal. Another example is the superconductiag quantum 
interferometric device (SQUID), which utilizes the quantum effects of superconductivity to 
detect magnetic fields. Other examples include imagery intelligence (IMINT), signals 
intelligence (SIGINT), electronic intelligence (ELINT) devices. In addition, active signal 
processing has been applied to femto-lasers. While these hardware-based systems provide 

1 5 enhanced signals in some circumstances, such systems are expensive to build and operate 
and can only physically manipulate or generate very specific types of signals. 

3. Stochastic Resonance Techniques 

Another exemplary active system or technique that has been implemented in 
20 software is stochastic resonance wherein noise is apphed to an image or a system being 
analyzed to improve at output signal-to-noise ratio (SNR). With stochastic resonance, a 
detector is modeled in terms of a bistable nonlinear dynamic system, which is enhanced by 
applying additive random noise element and a periodic sinusoidal forcing function. 
Stochastic resonance occurs when the SNR passes through a maximum as the noise level is 
25 increased. By properly selecting the amount of noise being applied, events of interest can 
be enhanced so that signals below backgroxind can be detected. 

Traditionally, stochastic resonance is employed using a classical noise source such 
as a standard Gaussian noise pattern. Alternatively, however, stochastic resonance may 
employed using some form of quantum noise. In this regard, it has be^ shown that 
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irradiation of a bistable system with a weak resonance laser field induces long-lived 
coherent tunneling motion whose amplitude exhibits a quantum stochastic resonance 
(QSR) as a function of the two-level-system relaxation rate and the fi.eld strength (See 
Makarov DE. Makri N. "Stochastic resonance and nonhnear response in double-quantum- 
5 well stractures" Physical Review B. 52(4):R2257-R2260, 1995 Jul 15). Essentially, QSR 
is an enhancement of the response of a driven quantum mechanical system by quantum 
noise. 

In any case, heretofore, traditional stochastic resonance and QSR have only been 
applied to point dynamical systems, i.e. to systems providing only a signal data point 
1 0 varying in time. One other technique for enhancing stochsistic resonance is so-called array- 
enhanced stochastic resonance (AESR), which extends stochastic resonance to 2-D 
dynamical systems, though still limited to static systems. Althotigh an improvement, 
AESR is typically limited to the analysis of systems with relatively few degrees of 
fi-eedom. 

i 5 Thus, most, conventional active signal analysis systems are not implemented in 

software and those that have been implemented in software typically have significant 
limitations. Accordingly, it would be desirable provide new or improved techniques for 
implementing active signal processing in software and it is to these ends that aspects of 
invention are drawn. By implementing active signal analysis in software, all of the 

20 advantages associated with software systems over hardware systems can be exploited, 
including lower cost, increased flexibility, increased robustness, portability, etc. 



implemented in software capable of detecting events of interest in spatial, temporal or 
spatio-temporal systems within a combinatorial state space wherein the state space is 
25 discrete and continuous and wherein background noise is nonstationary. Ideally, the 
system would be capable of analyzing signals received from a wide variety of arrayed 
platform detector configurations including linear, arrayed, spectral, and time-varying and 
can handle iV-attribute systems. Moreover, ideally, the improved system would be capable 
of analyzing signals with a signal to background ratio of 10*^. Insofar as the analysis of 



More specifically, it is desirable to provide an active signal processing system 
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biological samples applied to a microarray platform is concerned, fhe syst^ or technique 
would ideally be capable of detecting and quantitating specific hybridization (e.g., genomic 
events) and specific binding events (e.g., proteomic events). Insofer as genomics is 
concerned, the system should be capable of processing data fi*om a high-density gene 
5 expression microarrays (i.e. array is providing 10^ to 10^ spots per array) and which 
represents discrete, spatially indexed systems having high probe densities (i.e. from 10^ to 
10^ probes per spot). 

Summary of the Invention 

10 In accordance with the invention, techniques are provided for performing active 

interferometric signal analysis in software. Depending upon the implementation, the 
techniques provide for both detection and qioantitation analysis by exploiting constructive 
or destructive interferometric analysis (or a combination of constructive and destructive 
interferometric analysis) using reverberant convergence to detect resonance events. The 

15 techniques achieve software emulation of wave-particle interactions and wave-wave 
interactions and can operate in either the fi-equency domain or the phase domain. The 
techniques may be used for analyzing static spatial systems, static data from arrayed 
measurement platforms, dynamical systems, spatio-temporal systems or plasma systems. 

Briefly, the techniques exploit expressor functions designed to reject any 

20 interfering noise or background clutter from possible events of interest within arrayed data 
to be analyzed. More specifically, the expressor functions are designed so as to extract 
spectral invariants of possible events of interest associated with the array device used to 
detect the arrayed data. Various aspects of the invention are directed to the generation of 
the expressor functions. Typically, the expressor functions are designed to incorporate 

25 some form of non-Gaussian noise, such as quantum noise, pseudorandom noise, 
polycyclostationary noise, cyclostationary noise, stationary noise, non-stationary noise or a 
systemic bias. Preferably, prior to the application of the expressor function to arrayed data, 
the array is pre-conditioned so as to convert the arrayed data to a spectral domain in which 
spectral harmonics parameterize events of interest to a pre-determined dynamical system. 
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Preconditioning may be performed by applying a preconditioning function incorporating, 
e.g., a 1-D Fourier function, a 2-D Fourier function, an N-D Fourier function, a time 
division multiplexing (TDM) function, a wavelength division multiplexing (WDM) 
function, a firequency division multiplexing (FDM) function, a radial bias function, a 
5 wavelet kernel fimction, a fractal function, or a soliton function. The arrayed data itself 
may be the form of a spatial 1-D array, a spatial 2-D array, a spatial point emitter array, a 
temporal point emitter array, a spatio-temporal point emitter array, a spectral point emitted 
array or a virtual array constructed by combining spatial separate point emitters. 

In accordance with an exemplary method, a technique is provided for actively 

10 analyzing a signal pattern representative of arrayed data to identify events of interest 
therein. A signal pattern representative of arrayed data is input and resonance patterns are 
generated based on interference between synthetic noise and the signal pattern. Next, 
resonances within the resonance patterns associated with events of interest are detected. 
Preferably, the signal pattern is preconditioned prior to the step of applying synthetic noise 

15 to the signal pattern. 

In a first exemplary system embodiment, also referred to herein as an "open loop 
interferometric system", a system is provided for analyzing an arrayed signal pattern 
generated by an arrayed platform device to identify events of interest within the signal 
pattern. The system includes an expressor function input unit , for inputting previously 

20 generated expressor fimctions capable of extracting spatial, spatio-temporal or spectral 
invariants of events of interest associated with the particular arrayed platform device being 
used. The system also includes a preconditioning unit for preconditioning the arrayed 
signal pattern. An active interferometric coupler operates to convolve the preconditioned 
arrayed signal pattern and the expressor functions so as to interferometrically enhance 

25 portions of the preconditioned pattern associated with events of interest, if any, present 
within preconditioned pattern. The system also includes a resonance marker detector for 
identifying the occurrence of events of interest within the convolved signal pattern. 

In a second exemplary system embodiment, also referred to herein as a "goal- 
directed interferometric system", the system for analyzing an arrayed signal pattem 
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includes an expressor function input unit for inputting previously-generated expressor 
functions and an adaptive interferometric coupler for preconditioning the arrayed signal 
pattern so as to convert the arrayed signal pattern to a spectral domain while 
simultaneously convolving the signal pattern and the expressor functions so as to 
5 interferometrically enhance portions of the signal pattern to identify events of interest, if 
any, within the signal pattern. 

In a third exemplary system embodiment, also referred to herein as a "self- 
organizing interferometric system", the system includes a preconditioner for 
preconditioning the arrayed signal pattern and an expressor function adaptation unit for 
10 generating preconditioned expressor functions based on input canonical expressor 
functions (i.e. generaV expressor ftmction not associated with the particular platform being 
used) and based on a preconditioned signal pattern while interferometrically enhancing 
portions the preconditioned signal pattern so as to identify events of interest, if any, within 
the signal pattern. 

15 In a fourth exemplary system embodiment, also referred to herein as an "iterative 

interferometric system", the system includes an expressor function input unit for inputting 
previously-generated expressor fiinctions derived for the platform being used, a 
preconditioner for preconditioning the arrayed signal pattern and a iterative convolution 
coupler. The system also includes an adaptive controller unit for controlling the coupler to 

20 iteratively and selectively convolve the expriessor functions to a preconditioned signal 
pattern until a predetermined degree of convergence is achieved so as to identify events of 
interest within the enhanced signal pattern. 

A variety of techniques are set forth for performing active signal analysis in- 
software. Method, system and computer product embodiments are provided. 



The features, objects, and advantages of the present invention will become more 
apparent fix)m the detailed description set forth below when taken in conjimction with the 
drawings in which like reference characters identify correspondingly throughout and 



25 



Brief Description of the Drawings 
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wherein: 

FIG. 1 is a block diagram providing an overview of system components of the first 
exemplary embodiment of the invention (the open loop interferometric system), which 
provides sq>arate expressor function generation imit, preconditioner, active interferometric 
5 coupler and resonant marker detector components; 

FIG. 2 is a flowchart providing an overview of method steps for the first exemplary 
embodiment of the invention; 

FIG. 3 is a block diagram providing an overview of system components of the 
second exemplary embodiment of the invention (the goal-directed interferometric system), 
1 0 which provides an adaptive interferometric coupler; 

FIG. 4 is a flowchart providing an overview of method steps for the second 
exemplary embodiment of the invention; 

FIG. 5 is a block diagram providing an overview of system components of the third 
exemplary embodiment of the invention (the self-organizing mterferometric system), 
1 5 which provides a expressor fimction adaptation unit; 

FIG. 6 is a flowchart providing an overview of method steps for the third 
exemplary embodiment of the invention; 

FIG. 7 is a block diagram providing an overview of system components of the 
fourth exemplary embodiment of the invention (the iterative interferometric system), 
20 which provides a adaptive controller for iteratively controlling a separate iterative 
interferometric coupler; 

FIG. 8 is a flowchart providing an overview of method steps for the fourth 
exemplary embodiment of the invention; 

FIG. 9 is a graph illustmting combinations of arrayed data signal patterns that may 
25 be processed by any of the embodiments of the invention; and 

FIG. 10 is a graph illustrating combinations of preconditioning functions and 
expressor functions that may be used in coimection with any of the embodiments of the 
invention; 

FIG. 11 is a graph illustrating an exemplary platform array detector and active 
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interferometric analysis system that may e^loit the invention; 

FIG. 12 is a gr^h illustrating various platform array detectors that may be used in 
connection with any of the embodiments of the invention; 

FIG. 13 is a gr^h illustrating the structure of an exemplary microarray platform 
that may be used to generate arrayed data for use by any of the embodiments of the 
invention; and 

FIG. 14 is a flowchart stmmiarizing the overall method of the invention, at least in 
accordance with an exemplary embodiment thereof. 
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Detailed Description of Exemplary Embodiments 
With reference to the remaining figures, exemplary embodiments of the invention 
will now be described. 

5 A. OPEN LOOP INTERFEROMETMC SYSTEM 
1. Overview 

FIG. 1 illustrates, at a high-level, exemplary system components for the first 
exemplary implementation of the invention. Components employed only during a design 
phase are shown on the left; components employed during the actual analysis of an input 

10 signal are shown on the right. Briefly, during the design phase, expressor function 
generation unit 100 generates one or more designer expressor functions based upon the 
specific characteristics of a platform array detector 102 to be used to detect signal patterns 
to be processed and also based on calibration data derived from platform array 102. The 
calibration data is generated based on known events of interest That is, input signals or 

15 physical samples containing known events of interest are applied to the detector to 
generate the calibration data. In one example, the platform array detector is a genomic 
biochip/microarray and the calibration data is derived from a biological sample containing 
known gene expressions. The characteristics of the platform array detector that are input 
by the expressor fimction generation unit include both the layout of the array as well as its 

20 coherent noise or background signal characteristics. Once the expressor functions have 
been generated for the particular platform array detector based on the calibration data, the 
design phase is complete and, thereafter, the expressor function generation unit is no 
longer required unless an alternative platform array detector is employed. 

Components employed during the actual operational phase of the system include 

25 the aforementioned platform array detector 102 that generates a detected signal pattem and 
a preconditioner unit 104 that preconditions the detected signal pattem so as to convert the 
signal pattem to a spectral domain (for an example of a technique for transforming the 
signal pattem to a spectral domain, please see U.S. Pat App. Ser. No. 10/430,664 entitled 
"Method and System for Characterizing Microarray Output Data" fiLled on May 5, 2003 
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which is hereby folly incorporated by reference) in which spectral harmonics parameterize 
events of interest to a predetermined dynamical system (for which the applicable 
characteristics are input into the preconditioner unit The operational phase components 
also include an active interferometric coupler 106 that convolves the preconditioned signal 
5 pattern and the previously generated expressor functions. The convolution is performed so 
as to interferometrically enhance portions of the preconditioned signal pattern associated 
with events of interest, if any, that are present within the preconditioned signal pattern. As 
explained more folly below, the convolution is performed via reverberant convergence so 
as to emulate active interferometric enhancement. Finally, the system also includes a 

1 0 resonant marker detector 108 that processes the convolved signal pattern so as to identify 
particular events of interest, if any, appearing within the convolved signal pattern. 

With the exception of the platform array detector, which is a physical hardware 
component, all other- components illustrated in FIG. 1 may be implemented in software, 
hardware, firmware or some combination thereof. In particular, each component may be a 

15 computer product operative to perform the functions described. In one example, each 
component is a software module operating within a single generally programmable 
computer. In other cases, the software modules are implemented using separate 
microprocessors or application specific integrated chips (ASICs). In still other examples, 
some of the components are implanted in software whereas others are completely hard- 

20 wired. As can be appreciated, a wide range of possible implementations is consistent with 
the invention and no attempt is made herein to describe all possible configurations. 

Method steps associated with the open loop interferometric system of FIG. 1 will 
now briefly be simmiarized with reference to FIG. 2. Although these steps are 
advantageously performed by the system of FIG. 1, the steps can be performed by any 

25 other suitable system. Conversely, the sj^tem of FIG. 1, or portions thereof, may be 
employed to perform methods different fi-om that which is illustrated in FIG. 2. The same 
is true for the various otiier figures of the patent application. In any case, at step 110, a 
signal pattem derived firom a particular platform array detector is input and, at step 112, 
expressor fonctions designed to extract spectral invariants for events of interest detectable 



wo 2004/102456 PCTAJS2003/021525 

-13- 

by the particular platform array are input Then, at step 114, the input signal pattern is 
preconditioned to convert lie signal pattern to the spectral domain in which spectral 
harmonics parameterize events of interest to the aforementioned predetermined dynamical 
system. At step 1 16, the preconditioned signal pattem and the expressor functions are then 
5 convolved using reverberant convergence so as to emulate an active interferometric 
enhancement of portions of the signal pattern associated with events of mterest thus 
yielding a convolved signal pattem (also referred to herein as an enhanced signal pattem. 
Fmally, at step 1 18, the convolved signal pattem is examined to identify events of interest 
therein. 

10 It should be noted that a hst of specific events of interest need not be input at step 

1 12. Esther, the technique operates to detect possible events of interest present in the input 
pattem based solely on the input pattem and the previously-generated expressor fiinctions. 
Hence, the technique of FIG. 2 is not merely a diagnostic technique that seeks to 
determine whether particular pre-determined events of interest are present in the input 

15 signal pattem. This will be more readily apparent from the detailed examples which 
follow. 

2. Exemplary Implementation 

The mathematical equations performed by the various components of FIG. 1 or 
20 during the various method steps of FIG. 2 will now the set forth in detail. In this example, 
quantum expressor functions are employed. However, as described below, otiier expressor 
functions can alternatively be used, with" the equations modified as needed dependent upon 
factors such as signal to background strength, extraction core length of the events of 
^ interest bemg analyzed, platform coherence, and desired sensitivity. 
25 The exemplary equations are particularly relevant for conditions where the 

background is 10 to 1000 times stronger than the signal. 

Initially, during the design phase, expressor function genemtion vmit 100 of FIG. 1 

inputs platform array detector characteristics represented mathematically as: 



wo 2004/102456 



-14- 



PCTAJS2003/021525 



Y[a(i\^,Af,r,A)] 

where Q,{NyM,T,A) denotes a preconditioned extraction core corresponding to 
one or more events of interest The indices N (0:n), M(0:m), T(l:k), A(0:p) respectively 
refer to the spatial, temporal and spectral dimensions of the physical or virtual arrayed 
platform described in FIG 11. Details of preconditioning the extraction core corresponding 
to event of interest are given in U.S. Pat App. Ser. No. 10/430,664 entitled "Method and 
System for Characterizing Microarray Output Data" filed on May 5, 2003. 

The Expressor Function Generation unit also inputs calibrated data (j exemplars) 
represented mathematically as: 

J 

Based on the input calibrated data, the expressor function generation unit performs 
calculations provided below. 

In some embodiments. Quantum Expressor Functions (QEF*s) are calculated based 

upon the signal pattern based upon one or more events of interest The QEF is generated 

by the steps of: 



(a) 


selecting a canonical coupled Spin Boson basis system by first 




conditioning a Hamiltonian for that system; 


(b) 


calculating harmonic amplitudes |Pm| fot the Hamiltonian 


(c) 


generating an order fimction (OF); 


(d) 


measuring the entrainment states of the OF using exemplars with groxmd 




tmth data; and 


(e) 


modulating the OF using calibrated exemplars to jrteld the QEF 
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Fully characterized exemplars of the event of interest are used to develop a 
platform specific QEF. Hamlltonians are utilized as fhey provide a mathematical basis for 
determining the dynamical behavior of a system at each coordinate and at each momentum. 
5 Hamiltonian systems are fiirther characterized by kinetic and potential energy. A key 
property of Hamiltonian systems is that they are conservative and have no dissipation of 
energy during relaxation. The use of a Spin Boson Hamiltonian, which is a particular type 
of Hamiltonian systems, permits the exploitation of quantum stochastic resonance 
phenomenonology. Hamilton's Equations and Hamiltonian functions pertain to well 

1 0 understood concepts in dynamical systems theory and dynamical equations of motion. (See 
"Tlieory and Problems of Lagrangian Dynamics", Ed. Dare A. Wells, Schaum's Outiine 
Series, McGraw-Hill Company, NY, 1967, pages 316-321; Elements of Hamiltonian 
Mechanics", Dter Haar, Pergamnon Press, 1960, pages 95-102; and "Statistical 
Mechanics", R.K, Pathria, Pergamon Press, Toronto, 1972, pages 6-7, 136-147). 

15 The QEF encodes quantum noise to drive order in an asymmetric bistable 

quantum-mechanical system. The appearance of a resonance requires an asymmetiy in the 
energies of the two states. A rate equation constructed for transitioning between the two 
states of a bistable system, such that the dynamics can be characterized in terms of 
transition rates O + and O- between two asymmetric quantum superposition states where 

20 the drive fi:equency and interwell transition rates are much slower than intrawell relaxation 
rates. The signal to noise ratio (SNR) of such a superposition system is given by: 
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SNR = 



4 l+exp[^o/^ii^o] 




Where Kb is the Boltzmaxm's constant, T is the temperature, and the sinusoidally 
modulated asymmetry energy e is given by: 



For the above bistable system, the power spectrum S(a)) may be represented by the Fourier 
transform of C(x), containing a roughly Lorentzian broadband noise background and 5- 
function peaks at <d = 0, the driving frequency (Os, and its harmonics. A measured 
correlation function for the quantum noise is given by C(t) = <n+i (t)n+i (t¥ t)>, where each 
10 n+i = 0 or 1, and the <, indicates an average over t over many data points i taken at equal 
. intervals^ and is given by 

The probability of being in the + quantum state at / after being in the state go is 
given by 



C(T) = (i 



5 



oo). 



15 



A custom Hamiltonian, couples the above system to an ensemble of harmonic 



oscillators is given by: 
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In the above Hamiltonian, e denotes the asymmetric energy, A is the tunneling 
matrix element, and Oi are Pauli spin matrices, and Ct) is a harmonic oscillator creation 
5 operator witti frequencies cOr,. 

An important aspect is to couple the transformed and preconditioned discrete 
microarray output to a mathematical model for a quantum-mechanical dynamical system 
with specific properties. Specific exemplary parameters for use in calculating the 
Hamiltonian are those proposed by A J. Legett et al., Kevins of Modem Physics^ 59, 1, 
10 1987 and A.O. Caldiera and A J. LcgeU Annals of Physics, 149, 374, 1983. The parameters 
are important for an offline simulation of this Spin Boson system on a digital computer. 
The empirical observables are then collected and used to estimate and compute spectral 
properties, which are actually used by the method. 

Experimental or analytical vaUd parameters for fte above system will work with 
15 the technique because the robustness of the method depends only on the bulk and 
qualitative properties shown by this system and the properties of its power spectral 
harmonics. This is an important because the method is substantially immune to nuances 
and specifies of the actual driving ma&ematical system. 

The information about the effects of the environment is contained in the spectral 
20 density J(aJ^ 
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The hannonic amplitudes determine the weights of 5 spikes of an av^ged spectral 
power density in an asymptotic state S"((i)). e refers to the coupling strength between the 
spin and boson systems, and Po defines the equilibrium state in the absence of driving 
5 force. Power Spectral Density (PSD) is a well understood concept is signal processing and 
basis for many computational algorithms. For a spectral signature, partitioned into bands of 
interest, power spectral density is formally as the set of measurements of average power in 
each spectral band, normalized by bandwiddis. Formal definition can be foimd in 
"Cyclostationarity in Communications and Signal Processing", Ed, William Gardner, 
1 0 1994, ffiEE Press, NY, pages 46-47. 

The power amplitudes rj^ in the /nth jfrequency component of asymptotic state 

space are calculated using 

7„(£2,^) = 4;r|P,(Q,f)| 



15 



and the phase shift is given by 



:^^(Q,e ) = tan- 



lm(p^(Q,g)r 
Re/>^(£2,^) 



The analytic for the external force is given by 
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The parameters y, eo are predetermined and are design specific. Typically, values of 
0.001 and 0.0001 are used for eorcspectively. In the above e^qptession \P„\ determine 

5 the weights of the S spikes of the averaged spectral power density. 

For some embodiments, the QEF is designed by matching the power spectral 
density (PSD) amplitude of preconditioned, extracted events of interest to that of the Spin 
Boson system described above so that stochastic and deterministic time scales match and 
so that the time scales couple back to noise statistics and degree of asymmetry. The 
10 technique employs a fully automated iterative conjugate gradient relaxation method for 
spectral matching between asymmetric base system and the transformed, preconditioned 
extraction core. The determination of the QEF depends on the specifics of bioelectronics 
substrates used for actual analysis. The method is however generalizable to all arrayed 
embodiments. In addition, the melliod is highly scalable to array dimensions (as the 
15 offline design trade-space time does not matter to computational complexity). As the 
system is an overdetermined coupled system, convergence criteria and stability of 
relaxation metiiod directly relates to downstream resonance effectiveness. 
Generation of Order Function 

The order function (OF) of ground truth is generated as follows. The order 
20 fimction (OF) is for groimd truth wherein ground truth represents a state wherein a positive 
signal to noise ratio (SNR) is expected for pixel intensities of selected exemplars used for 
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developing an QEF. 

A classical Order Function (OF) is used with the following form: 

1 • 

where i= -/^ and \^\deh{0)e'^^ and 0<^<1. 



0 



The term 

hf^ ^jd6hi0)e'^'^^ is referred to as the Diado integral. 

0 

10 In this application, QEF are designed from a canonical dynamical system drawn 

from a family of exponential systems. Daido integral is a specific term used to construct an 
Order Function, which is then used to construct the QEF. 

The OF is calculated using an order match which impUes that the variance of the 
density function for a specific exponential family approaches 0 (within 0.000001 - 0.0001). 

1 5 The free energy for the density function is given by 



where and represent state vectors. Jf * represents the random observable in 
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symmetric bilinear form, and X denotes die characteristic function, calculated using 
extractions from exemplars with reference objects of interest known (see D.C. Brody and 
L.P. Hughston, "Geometry of Quantum Statistical Inference", Physics Review Letters, 
77(14), pp. 2851-2855, 1996). 
5 Determination of Entr ft^inniftnt States 

The OF derivation uses results from Daido's theory of multibranch entrainment of 
coupled nonlinear oscillators, wherein a number of different entrained states co-exist. For 
exemplars, an a priori measurement of entrainment is performed. This is done by using an 
approximate Daido IntegralQ with Zk = PSD maxima at regions where boimdary of event 
10 of interest > desired detection threshold (see Physics Review Letters, 77(7), 1406-1411.) 
Pre^multiplier constants are used to ensure that meets Ihe following tests: 

- maxima power spectrum density (PSD) matches; and 

- hi norm on even harmonics < e where e = NS-MRF Barrier assimiing complete 
spatial randomness (CSR) assumptions discussed in U.S. Pat. No. 6,136,541, entitled 

15 "Method and Apparatus for Analyzing Hybridized Biochip Patterns using Resonance 

Interactions Employing Quantum Expressor Functions all of which is hereby folly 

incorporated by reference**. Note that NS-MRF refers to the nonstationary Markov random 
« 

field representation. 

The notion of entrainment states is exploited, in part, because the method treats the 
20 spectral system embodied in an extraction core corresponding to the event of interest, as a 
special case of coupled nonlinear oscillators in equilibrium. However, due to device 
imperfections, signal degradation and other limitations, number of entrained binding states 
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(i.e., background and clutter conditions) may coexist. A PSD match is desired where the 
OF is Z-peaked (but single peaked around each event of interest centroid). Absence of a 
single peak implies imperfect or lack of presence of an event of interest It also defines the 
resonance loci for Ihis method (i.e., where maximum SNR enhancement) is obtained. 



The OF of ground truth is modulated to yield the QEF as follows. Under controlled 
calibration, as stated above, maximal SNR enhancement (optimal resonance) is achieved 
when OF yields a single peak. It is an important design point for matching PSD of 
coupling Spin Boson system to the synthetic QEF. The specific form of the QEF to be 
10 used is the parameters of a generic OF shown above. So the exemplary method exploits 
two connotations of OF: (a) parametric form for the QEF (that is closer to the classical 
form) and (b) as exponential attractor for a dissipative system. The two OF's are then 
recoupled and convolved with spectrally transformed exemplar used for constructing the 
QEF. 

1 5 The resulting QEF is given by: 



Typically, only first three to twenty-five PSD peaks are considered for spectral matching. 

Preferably the QEF is represented digitally using a matrix or array having the same 
number of elements in the spectrally transformed and preconditioned, extraction core to be 
20 analyzed. In the exemplary embodiment this pertains to tiie preprocessing done to the post- 
hybridization microarray intensity output to get it into a form where it can be convolved 
with the Quantum E?qpressor Function. 



5 



Ground Truth based Modulation of OF 
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A major limiting restriction in QSR that is avoided by the exemplary method 
pertams to matching the stochastic and deterministic time scales in "domain system" and 
the external coupling asymmetric dynamical system, since this has limited q>plicability to 
5 continuous data. By replacing the time scale match requirement witti ensemble spatial 
statistics (statistical mechanics of the signal emitters responsible for generating the pattern 
of interest) expressed via the aforementioned generators an entirely new class of analysis 
for discrete systems is enabled. By deriving the condition that signal emitters need to 
satisfy to maximize QEF generator capacity, a properly constructed and conditioned QEF 
10 will be able to extract and enhance an entire class of interest to produce one or more 
expressor functions represented mathematically as 

b. Preconditioning 

15 - 

The signal pattern output by platform array detector 102 is represented 

matiiematically as: 

20 
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Using tiie input signal pattern, the Preconditioner unit 104 performs the 
mathematical calculations detailed in U.S. Pat App. Ser. No. 10/430,664 entitled "Method 
and System for Characterizing Microarray Output Data" filed on May 5, 2003 which is 
hereby incorporated by reference, to generate the preconditioned output signal pattern 
5 represented mathematically as: 



10 



c. Convolution Coupling 
Active interferometric coupler 106 performs the following mathematical 
calculations using the preconditioned signal pattern: 



-7:(0) 




15 



is defined as a vector containing the preconditioned components from an event 
of interest. 




(0 




where represents is the QEF after / convolutions. 



20 




(0)— (1) 




-t(1)t::(2) 



^(0)— (2)— (1) 




where 



represents 
quantity) 



the QEF developed in the preceding step and (its dimensionless 
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and where !b5 represents die i-tfa pertuibation to the QEF, induced by perturbing one 
of its spectral components. 



for k = 1 to w 

for j= 1 to 1000 (set to a large counter value) 

perturb the component of QEF as below 



(k) = [&-\k) + jC, sxixiwj + CJt 



where 



r X if X>0^ 



and 



C =if^) 

1 3 3 60 * ^ denotes a small constant; 

W = 

Let 0 the variance computed from the values 



/ 



pc 



][^CIX(^^ ) where fpc denotes the preconditioned spectral vector 

corresponding to a known event of interest present in the arrayed pattern 
being analyzed. 

2 

As an example, fpc refers to the spectral components of the positive control. 
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The convolution iteration can be e^qpressed as: 



10 



25 



where fnc refers to the spectral components of a canonical negative control, or 
preconditioned footprint of an event of interest known to be absent in the arrayed 
image. 

After each convolution iteration check for monotonicity of 



if yes — > then exit loop to perform global QEF iterations (that means this 

particular k component is important, i.e we are diverging from 
15 the negative control.) 

If no — > dien continue 
end J loop 
end k loop 

20 Global QEF iterations may be performed if monotonic divergence is detected 

between tihe preconditioned extraction core being analyzed and the canonical negative 
control, then the same convolution coupling operations are r^eated for all the spectral 
harmonics. The global QEF iterations are provided by: 



Set VZ 

For ^ — 1 to 25 (chosen to be a small count value) 
Compute 
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f Nc_ 

where ^1 c^tures the variance of Uie components of 2^Clx(^f ) 



and 

C2 = Ci+£ 



'^Parseval Avg from Pos. Con, PM^ 
Parseval Avg from Neg. Con PM y 



where Parseval Avg. firom Pos. Con. PM refers to the parseval number for a canonical 
event of interest known to be present, and Parseval Avg. from Neg. Con. PM refers to the 
parseval number for a canonical event of interest known to be absent 

10 

and ^ is chosen to be small, 0.0001 . 

Again after each coiipler iteration compute the term 



Successively compute 



After each convolution iteration check if 

^^^>l AND ^^>1 

If the conditions of the above test are met, resonance is concluded and event of interest is 
called present 
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If the monotonicity test fails, them the preconditioned test pattern is nonnalized using the 
expressions below. 



^y-O)_ya-0>_^(.)_y;(.) 



(J) 70) 



fhen 



V 



pc J Nc for any component 

25 J 



10 



The above assumes an object with 25 spectral harmonics of interest 



15 



The detailed equations for the coupler unit and resonance detector unit are given below: 



for 



y = 1 to N 
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10 



i = k 



A / 



i = k 



2 

^ y-i,/+2 



Check 



Ay + l->y + 2 



>1 



>1 



>1 



^ j 



^ 1 2;^ 
0 = 

3 360 



Where ^ = f i^l) 



and 



V 



^man PSD 



pm,njc y 
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^ = ^c, 0 • {sin(w, (r) + (p^, )} 



{=0.0001} 



PSD 



QEF pc.pm 



V 



max(PSDQEp.p,p^) ^ 



10 



PSD 



QEF,NC,pm 



max(PSDQEp,,e,pm) J 



15 



d. Resonant Marker Identification 

Finally, resonance marker detector 108 performs the following 
mathematical calculations using ifae convolved signal pattern to idoitify the events 
of interest within the convolved signal pattern. 



20 
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The resonant iteration is teiminated when 




>1 



a/ -^7 + 1 



>1 



A7^y + 1 



A7 j 



or when iteration counter / exceeds preset "JV" (e.g., 10^ iterations) (for digital 
approximation to analog dynamics). 

Preferably, the resonance interaction is performed digitally by applying a matrix 



stimulus in combination with a matrix representative of the preconditioned extraction of 
event of interest. 

The final result of interferometric analysis is the readout of locations wherein 
resonance has occurred and as identified by row and column number (/j). 



For example, with regard to gene expression, once resonances peak are observed by 
applying a QEF to specific object extractions, the object addresses ifcJC) of those locations 
are mapped into the physical feature or gene table (which contains the complimentary 
RNA or cDNA probe sequences associated witii object locations) to thereby identify the 
gene being e7q)ressed, if any, present in the sample being analyzed. This is a simple table 
look-up resulting in a direct readout of the names of the expressed genes. 



representative of the resonance equation to a matrix representative of the resonance 



Resonance Output Interoretation to Identify Gene Expression 
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To suxmnarize the open loop interferometric technique of FIGS. 1 and 2, the 
5 preconditioning of arrayed data to the spectral domain is performed independently prior to 
coupling with the expressor function. Hence, the reverberation convergence performed by 
the active interferometric coupler to achieve a resonance state is performed independent of 
preconditioning. The reverberations are an open loop process terminated upon detection of 
a predetennined condition (i.e. observation of a resonance marker.) The active 
10 interferometric coupler implements either a destructive or a constructive interference of 
two dynamical systems. 

B. GOAL-DIRECTED INTERFEROMETRIC SYSTEM 

FIG. 3 illustrates exemplary system components for the second exemplary 

15 implementation of die invention. As before, components employed only during a design 
phase are shown on the left and components employed during die actual analysis of an 
input signal are shown on the rigjit Some of the components of the system of FIG. 3 are 
the same as those shown in FIG. 1 and hence will be not described again in detail Again, 
during a design phase, an ejqpressor function generation unit 200 generates expressor 

20 functions based upon calibration data and based on the specific characteristics of a 
platform anay detector 202 and during the operational phase, the platform array detector 
generates a signal pattern for analysis. However, no separate preconditioner imit is 
provided. Nor is any separate resonant maker detector provided. Rather, a single adaptive 
interferometric coupler 204 is provided that operates to simultaneously precondition the 

25 detected signal pattern to convert the arrayed signal pattem to a spectral domain while 
convolving the signal pattem and the expressor fiinctions to interferometrically enhance of 
portions of the signal pattem and to identify events of interest, if any, within the 
preconditioned signal pattem. Again, with the exception of the platform array detector, 
which is a physical hardware component, all other components illustrated in FIG. 3 may 
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be implemented in software, hardware, firmware or some combination thereof. 

Note that adaptive interferometric coupler 204 does not merely perform the 
operations previously described with respect components 104, 106 and 108 of FIG. 1. 
Rather, adaptive interferometric coupler 204 perfiDrms a different set of operations. 
5 Briefly, with the technique of FIG. 3, coupling intelligence is merged with preconditioning 
and the expressor functions are pre-coupled to the preconditioned functions, as will be 
apparent from the detailed example below. Convergent reverberation is performed 
simultaueously with preconditioning, which is now a closed loop process. Hence, with the 
technique of FIG. 3, preconditioning does not merely transform spatial arrayed data to 
10 spectral data but also transforms the data to a pre-determined state of a predetermined 
dynamical system. If that state is reached, then resonaoce is implicitly assumed to have 
occurred. Thus, completion of the adaptive interferometric coupling procedxire implies the 
detection of resonant markers and so separate resonant marker detection unit is not 
required. 

15 Method steps associated with the goal-directed interferometric system will now 

briefly be simimarized with reference to FIG. 4. At step 206, a signal pattern derived from 
a particular platform array detector is input and, at step 208, expressor functions designed 
to extract spectral invariants for events of interest detectable by the particular platform 
array detector are input. Then, at step 210, the input signal pattern is preconditioned to 

20 convert the signal pattern to the spectral domain while convolving ttie signal with the 
expressor ftmctions so as to emulate an active interferometric enhancement of portions of 
the signal pattern so as to identify events of interest therein. As with the technique 
described above, a hst of specific events of interest need not be input during the 
operational phase. Rather, the technique operates to detect possible events of interest 

25 present in the input pattern based solely on the input pattern and the previously-generated 
e>q>ressor functions. 

One of ordinary skill in the art will appreciate that the equations provided above in 
coimection with the open loop interferometric system may be accordingly adjusted such 
that the resonant marker unit and the convolution coupling unit may be combined so that 
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the threshold check for resonance is performed during the post-convolution normalization. 



wo 2004/102456 PCTAJS2003/021525 

-35- 



C, SELF-ORGANIZING INTERPEROMETRIC SYSTEM 

VIG. 5 illustrates exemplary system components for the third exemplary 
implementation of the invention. Unlike the systems thus far described, no separate design 
5 phase is required. Rather, all components operate during the actual analysis of an input 
signal. Ag^in, some components are the same as those described above and hence will be 
not described again in detail. In use, a platform array detector 302 generates a signal 
pattern for analysis. The signal pattem is preconditioned by a preconditioner unit 304, 
which, as before, operates to convert the signal pattern to a spectral domain in which 

10 spectral harmonics parameterize events of interest to a predetermined dynamical system. 
The preconditioned signal pattem as well as characteristics of the platform array detector, 
calibration data and a set of canonical expressor functions are then all input by an 
expressor function adaptation unit 306. The canonical expressor functions are generalized 
expressor functions that have not been tailored to the particular platform being used. No 

15 separate design phase expressor fimction generation unit is provided. Nor is a separate 
resonant maker detector provided. Rather, expressor function adaptation unit 306 operates 
to simultaneoiasly generate preconditioned expressor fimctions based on the preconditioned 
signal pattem and based on the canonical expressor functions while convolving the 
preconditioned signal pattem and the preconditioned expressor functions to 

20 interferometrically enhance portions of the signal pattem to identify events of interest As 
before, with the exception of the platform array detector, which is a physical hardware 
component, all odier components illustrated in FIG. 5 may be implemented in software, 
hardware, firmware or some combination thereof. 

Note that expressor function adaptation unit 306 does not merely perform the 

25 operations previously described with respect to components 100, 106 and 108 of FIG. 1. 
Rather, with the technique of the self-organizing interferometric system, the coupling 
intelligence and the reveiberation dynamics are pre-built into the expressor functions and 
so die expressor fimctions are not merely a representation of a spectral marker as with the 
technique of FIGS. I and 2 and with the technique of FIGS. 3 and 4. This will be more 
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readily apparent from the detailed example below. Also, although the preconditioning unit 
may be the same as in the technique of FIGS. 1 and 2, the preconditioned array now 
serves as a preconditioner to the expressor fimction to thereby produce a properly 
preconditioned e7q)Fessor, which itself corresponds to the detection of resonant marker. 
5 Hence, generating and outputting the properly preconditioned expressor function results in 
detection of the event of interest As the expressor functions are tailored to the 
preconditioned array, they cannot be used again to process other signal patterns detected by 
the platform array detector. Also, as will be described in greater detail below, with the 
technique of the self-organizing interferometric system, the operations of the expressor 

1 0 fimction adaptation unit operate in an internal closed loop. 

Method steps associated with the self-organizing interferometric system will now 
briefly be summarized witii reference to FIG. 6. At step 308, a signal pattem derived from 
a particular platform array detector platform is input and, at step 310, the signal patterns 
are preconditioned to convert the signal pattem to the spectral domain in which spectral 

15 harmonics parameterize events of interest to the eiforementioned predetermined dynamical 
system. At step 312, preconditioned expressor fimctions are then generated based on the 
preconditioned signal pattem and on the characteristics of the platform array detector and 
based on the events of interest This is performed using a closed loop process employing 
reverberant convergence so as to emulate active interferometric enhancement of portions 

20 of the original signal pattem associated with the events of interest. Hence, generation of 
the preconditioned expressor functions allows immediate detection of the events of interest 
without further processing. As before, no list of specific events of interest need be input 
during the operational phase. 



25 connection with the open loop interferometric system may be accordingly adjusted such 
that the modulation of the OF to generate a QEF is immediately succeeded by the 
convolution coupling. 



One of ordinary skill in the art will ^preciate that the equations provided above in - 
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D. ITERATIVE INTERPEROMETRIC SYSTEM 

FIG. 7 illustrates exemplary system components for the fourth exemplary 
implementation of the invention. As before, components employed only during a design 
phase are shown on the left and components employed during the actual analysis of an 
5 input signal are shown on the right and components akeady described will be not described 
Bjg^in in detail During the design phase, an expressor function generation unit 400 
generates expressor functions based upon calibration data and based on the specific 
characteristics of a platform array detector 402 and, during the operational phase, the 
platform array detector generates a signal pattern for analysis. The signal pattern is 

10 preconditioned by a preconditioner unit 404, which, as before, operates to convert the 
signal pattern to a spectral domain in which spectral harmonics parameterize events of 
interest to a predetermined dynamical system.' An iterative interferometric coupler 406 
convolves the preconditioned signal pattern and the previously-generated expressor 
functions so as to interferometrically enhance portions of the preconditioned signal pattern 

15 associated with events of interest that are present within the preconditioned signal pattern. 
However, unlike the technique of the open loop interferometric system of FIGS* 1 and 2, 
wherein the convolution coupler is a stand-alone unit that operates to generate a final 
convolved signal pattern, with the technique of FIG. 7, the iterative interferometric coupler 
operates under the control of an adaptive controller 408 to sequentially produce an 

20 iteratively convolved signal pattern. As explained more fully below, the adaptive 
controller operates to control the interferometric coupler to iteratively and selectively 
convolve expressor functions to the preconditioned signal pattern until a predetermined 
degree of convergence is achieved so as to allow identification of events of interest within 
the enhanced signal pattern. Again, with the exception of the platform array detector, all 

25 other components illustrated in FIG. 7 may be implemented in software, hardware, 
firmware or some combination thereof. 

Method steps associated with the iterative interferometric system will now briefly 
be summarized with reference to FIG. 8. At step 410, a signal pattern derived from a 
particular platform array detector is input and, at step 412, expressor functions designed to 
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extract spectral invariants for events of interest detectable by the particular platform array 
detector are input Then, at st^ 414, the input signal pattern is preconditioned to convert 
the signal pattern to the spectral domain spectral domain in which spectral harmonics 
parameterize events of interest to the aforementioned predetermined dynamical system. At 
5 step 416, the preconditioned signal pattern and the ejq>res5or functions are then iteratively 
convolved under the control of the controller to achieve a reverberant convergence so as to 
emulate active interferometric enhancement of the signal pattern associated and to allow 
detection of events of interest As before, no list of specific events of interest need be 
input during the operational phase. 
1 0 One of ordinary skill in the art will appreciate that the equations provided above in 

connection with the open loop interferometric system may be accordingly adjusted such 
that multiple QEFs may be employed to achieve a desired resonance level and the most 
relevant QEF is chosen after each convolution if a desired resonance is not achieved. 

15 E. GENERAL FEATURES 

Exemplary embodiments are described above wherein arrayed data generated by a 
platform array detector is processed to identify events of interest therein. A wide variety of 
types of arrayed data maybe processed. FIG. 9 lists specific examples of types of arrayed 
data that may be processed. The examples shown include spatial 2-D data, spatial 1-D 

20 data, point emitter data, temporal point emitter data, spatio-temporal point emitter data and 
spectral point emitter data. The different types of arrayed data may be separately processed 
or, as shown, combined to yield virtual arrays for processing. In FIG. 9, the symbol \jf 
generally represents a virtual array combiner operation, wherein individual spatial values 
are collected together into a single function or data array. For example, if sixteen 

25 individual point emitters are provided, the data from the sixteen point emitters is combined 
or merged to yield a single 4X4 virtual array (represented, e.g., as f(x,y)). The symbol i 
generally represents a temporal accumulation operation, wherein individual values are 
accumulated over time to yield a single time-varying function or data array. For example, 
if a siQgle point emitters is provided, the point emitter may be sampled periodically to 
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accumulate values to populate a time-varying function (represented, e.g., as fl[t)). In some 
cases, both virtual array and accumulation operations are performed to generate spatio- 
temporal data (represented, e.g., as f(x,y,t)). 
More specilBcally, FIG. 9 illustrates: 
6 (a) spatial data fix)m a 2«D spatial array 500 represented in tiie form of a 2-D 

function h(x,y); 

(b) spatial data from a 1-D spatial array 502 represented the form of a 1-D function 

k(x); 

(c) spatial data from a group of point emitters 504 (with each point emitteri, for i = 
10 1 ... n, providing a single non-time- varying value) assembled into a single 2-D function 

h(x,y) via a virtual array combiner 506; 

(d) temporal data from a single point emitter 508 accumulated over time via a 
temporal accimiulator 510 to yield a single 1-D function g(t); 

(e) spatio-temporal data from a group of individual point emitter arrays 512 (with 
15 each point emitter, for i = 1 ... n, providing a 2-D array non-time var5dng values) each 

s^arately accumulated over time via temporal accumulators 514 to yield individual 
functions p(n,y,t), which are in turn assembled into a 2-D spatio-temporal function via a 
virtual array combiner 516; and 

(f) spectral data from a group of individual spectral point emitters 518 (with each 
20 spectral point emitter, for i = 1 ... n, providing an output function f(X) varying with 

wavelength) assembled into a 3-D virtual spectral array via virtual array combiner 520, 
with values then accumulated over time via temporal accumulators 522 to yield a 4-D 
function g(x,y,A,,t). 

Additionally, the various types of arrayed data ((a) - (f)) may be further combined 
25 into a larger virtual array (g) by virtual array combiner 524. Also, any of the types of 
arrayed data ((a) - (g)) may be further filtered via linear or non-linear filters 526, to yield 
filtered array functions (h). The final collection of data is represented in the figure as 
arrayed data 528. As can be appreciated, a wide variety of arrayed data may be input, 
combined and manipulated to yield individual multi-dimensional input functions for 
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processing by the above-described systems. 

Also, a wide variety of preconditioner functions may be employed for use in 
preconditioning the arrayed data and a wide variety of expressor functions may be 
generated for use in processing the preconditioned arrays. FIG. 10 lists specific examples 
5 of types of expressor functions and types of preconditioning functions. The types of 
expressor functions are identified in terms of the type of noise associated with the 
expressor function. Examples of expressor functions include quantum expressor functions 
(QEFs), classical expressor functions (CEFs), classical statistical noise expressor functions 
(SCNEFs), pseudorandom noise expressor functions (PRNEFs), systemic bias expressor 

10 functions and expressor fimctions based on polycyclostationary noise (PCS), 
cyclostationary noise (CS), stationary noise (SN), non-stationary noise (NS), Examples of 
preconditioner functions include 1-D Fourier functions, 2-D Fourier functions, N-D 
Fourier functions, a time division multiplexing (TDM) functions, wavelength division 
multiplexing (WDM) functions, fi*equency division multiplexing (FD) functions, radial 

15 basis functions, wavelet kernel functions, fractal fimctions and soliton functions. In 
general, any combination of expressor function type and preconditioner function type 
illustrated in the figure may be employed. 

FIG. 11 illustrates an exemplary system 700 incorporating a biological platform 
array 702 that provides arrayed data to a software-based active interferometric analysis 

20 system 704. The platform array includes a soUd support 706 upon which various probes 
708 are attached via attachment chemistry portion 710. In the example of FIG. 11, the 
probes include one or more of cDNA, mRNA or oligonucleotide probes or proteins, 
peptides or oligopeptides. Following application of a biological sample, such as blood or 
urine, to the probes, a detector 712 detects a resulting pattern, which is then forwarded to 

25 the analysis system in the form of arrayed data. The analysis system exploits one of the 
general techniques described above to identify events of interest present in the biological 
sample. In this example, the events of interest are typically output in the form of one or 
more ''gene calls". 

A wide variety of different platform array detectors may be used to generate the 
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airayed data Examples of platforms detectors are listed in FIG. 12. The examples include 
optical platforms^ biomolecular platforms, ionic platforms, biomechanical platforms, 
optoelectronic platforms, radio frequency platforms, and other electronic microdevices. 
Exemplary biomolecular spatial array platforms include: hybridized spotted cDNA 
5 microarrays, synthesized oligonucleotide arrays, spotted oligonucleotide arrays, peptide 
nucleotide assays, single nucleotide polymorphism (SNP) arrays, carbohydrate arrays, 
glycoprotein arrays,, protein arrays, proteomic arrays, tissue arrays, antibody arrays, 
antigen arrays, bioassays, sequencing microarrays, sequencing by hybridization (SBH) 
microarrays, siRNA duplexes, RNAi arrays glass-based arrays, nylon membrane arrays, 

10 thin film arrays, polymer-substrate arrays, capillary electrophoresis arrays, genospectral 
arrays, electronic arrays, bead arrays, quantum dot arrays, glycan arrays, spotted wells, and 
six)tted well plates. Typically, any of the foregoing platforms may be implemented as a 
microarray. Components of a typical- microarray are shown in FIG, 13, which include a 
substrate, a probe portion, a label portion, a bioassay portion, and a detector. Exemplary 

1 5 substrates include glass, nylon, thin film, and polymer. Exemplary probes CDNA probes, 
spotted oligonucleotide probes and synthetic oligonucleotide probes. Exemplary labels 1- 
dye labels and 2-dye labels. Exemplary bioassay portions include prehybridization buffers, 
hybridization buffers and wash buffers. Exemplary detectors include laser scanners, 
confocal microscopy detectors, charge coupled devices (CCDs) and electronic readouts. 

20 The substrate has a certain layout represented by feature size, various controls (positive, 
negative dye and alignment), and probe density. 

FIG* 14 provides a sununary of an exemplary technique of the invention. Briefly, 
the figure illustrates method for actively analyzing a signal pattern representative of 
arrayed data to identify events of interest therein. At step 700, a signal pattern 

25 representative of arrayed data is input and, at step 702, resonance patterns are generated 
based on interference between synthetic noise and the signal pattern, in accordance with 
techniques already described. Thereafter, at step 704, resonances are detected within the 
resonance patterns associated with events of interest The syntiietic noise used at step 702 
is preferably in the form of a quantum expressor function, a classical e?q>ressor function. 
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classical statistical noise, pseudorandom noise, a systemic bias or some combination 
thereof. The arrayed data may be, as already explained, in the form of a spatial 2-D array, 
a spatial 1-D array, an N-D array, a temporal point emitter array, a spatio-temporal point 
emitter array, spectral point emitter array or a virtual array constructed by combining 



The description of the exemplary embodiments is provided to enable any person 
skilled in the art to make or use the present invention. Various modifications to these 
embodiments will be readily apparent to those skilled in the art and the generic principles 
defined herein may be applied to other embodiments without the use of the inventive 
1 0 faculty. Thus, the invention is not intended to be limited to the embodiments shown herein 
but is to be accorded the widest scope consistent with the principles and novel features 
disclosed herein. 



5 spatial separate point emitters. 
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CLAIMS: 

What is claimed is: 

1. A system for performing active interferometric analysis 
5 using an e?qpressor function to detect the presence of an event of interest 

within an arrayed signal pattern via a computationally induced interference 
mechanism. 

2. The system of claim 1 implemented in software. 

3. The system of claim 1 wherein the arrayed signal pattern 
1 0 comprises either static data or dynamic data. 

5. The system of claim 1 wherein the expressor function is a 
quantum e7q;)ressor Amotion. 

6. The system of claim 1 wherein the arrayed signal pattem is 
in the form of one or more of a spatial 1-D array, a spatial 2-D array, an N- 

15 D array, a tempoml point emitter array, a spatio-temporal point emitter 
array, spectral point emitter array, or a virtual array constructed by 
combining spatial separate point emitters. 

7. The system of claim 1 wherein the armyed signal pattem is 
generated &om a device chosen from the group comprising: optical 

20 platforms, biomolecular platforms, ionic platforms, biomechanical 
platforms, optoelectronic platforms, radio frequency platforms, electronic 
microdevices. 

8. The system of claim 1 wherein the arrayed signal pattem is 
generated from a device chosen from the group comprising: hybridized 

25 spotted cDNA microarrays, synthesized oligonucleotide arrays, spotted 
oUgonucleotide arrays, peptide nucleotide assays, single nucleotide 
polymoiphism (SNP) arrays, carbohydrate arrays, glycoprotein arrays,, 
protein arrays, proteomic arrays, tissue arrays, antibody arrays, antigen 
arrays, bioassays, sequencing microarrays, sequencing by hybridization 
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(SBH) mioroarrays, siKNA duplexes, KNAi arrays glass-based arrays, 
nylon membrane arrays, thin film arrays, polymer-substrate arrays, 
capiUary electrophoresis arrays, genospectral arrays, electronic arrays, 
bead arrays, quantum dot arrays, glycan arrays, spotted wells, spotted well 
5 plates. 

9. The system of claim 1 wherein the expressor function is 
designed to reject any interfering noise and backgroimd clutter from events 
of interest 

10. The system of claim 9 wherein the expressor function is 
1 0 designed so as to extract spectral invariants of events of interest associated 

with an array platform device used to detect the arrayed signal pattem. 

11. The system of claim 1 wherein the expressor function is 
comprised of frequency domain sequences, time domain sequences or 
spectral sequences, phase sequences or numeric sequences. 

15 12. A system for performing active interferometric signal 

analysis of static data from an arrayed platform in software. 

13. The system of claim 12 wherein the active interferometric 
analysis includes detection and quantitation analysis. 

14. The system of claim 12 wherein the active interferometric 
20 analysis includes constructive interferometric analysis. 

15. The system of claim 12 wherein the active interferometric 
analysis includes destructive interferometric analysis. 

16. The system of claim 12 wherein the active interferometric 
analysis includes software emulation of wave-particle interactions. 

25 17. The system of claim 12 wherein the active interferometric 

analysis includes software emulation of wave-wave interactions. 

18. The system of claim 12 wherein the active interferometric 
analysis includes one or more of frequency domain, time domain, or phase 
domain analysis. 
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19. The system of claim 12 wherein the active interferometric 
analysis exploits iterative convergence to detect resonance events. 

20. A system for performing active interferometric analysis by 
using reverberant convergence to detect resonance events. 

5 21 . The system of claim 20 exploiting quantum interferometric 

analysis. 

22. The system of claim 21 wherein the quantum 
interferometric analysis is applied to one or more of a static spatial system, 
a static data from arrayed measurement platforms, d3mamical systems, 

1 0 spatio-temporal systems and plasma systems. 

23. The system of claim 20 wherein the active interferometric 
analysis exploits non-classical noise^. 

24. The system of claim 23 wherein the active interferometric 
analysis exploits one or more of non-Gaussian noise, ergodic noise, and 

15 quantum-mechanical noise. 

25. The system of claim 20 wherein the active interferometric 
analysis exploits resonant interferometric signal analysis using a digital 
representation of spectral pulses. 

26. The system of claim 20 wherein the active interferometric 
20 signal analysis eisploits one or more of coherent, incoherent, synchronized, 

or asynchronous oscillations. 

27. The system of claini 20 wherein the active interferometric 
analysis elicits time domain or frequency domain convolution operations. 

28. Hie system of claim 20 wherein fh& active interferometric 
25 analysis e7q>loits expressor frinctions to detect the presence of an event of 

interest via a computationally induced interference mechanism. 

29. A method for actively analyzing a signal pattern 
representative of arrayed data to identify events of interest therein, the 
method comprising the steps of: 
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inputting a signal pattern representative of arrayed data; 

generating resonance patterns based on interference between 
synthetic noise and the signal pattern; and 

detecting resonances within the resonance patterns associated with 
5 events of interest 

30. Hie method of claim 29 wherein the synthetic noise is in the 
form of one or more of a quantum expressor function, a classical expressor 
fiinction, classical statistical noise, pseudorandom noise, or a systentiic 
bias. 

10 31. Hie method of claim 29 wherein the arrayed data is in the 

form of one or more of a spatial 2-D array, a spatial 1-D array, an N-D 
array, a temporal point emitter array, a spatio-temporal point emitter array, 
spectral point emitter array or a virtual array constructed by combining 
spatial separate point emitters. 

15 32. The method of claim 29 further including the step of pre- 

conditioning the signal pattern prior to tibe step of generating resonance 
patterns. 

33. The method of claim 32 wherein the step of pre- 
conditioning the signal pattern is performed by applying one or more 

20 preconditioning functions in the form of a 1-D Fourier function, a 2-D 
Fourier fimction, an N-D Fourier function, a time division multiplexing 
(TDM) function, a wavelength division multiplexing (WDM) fimction , a 
frequency division multiplexing (FD) function, a radial basis fimction, a 
wavelet kernel function, a fractal function, a soliton fimction. 

25 34. The method of claim 29 wherein the arrayed data is 

generated from a device chosen from the group comprising: optical 
platforms, biomolecular platforms, ionic platforms, biomechanical 
platforms, optoelectronic platforms, radio frequency platforms, electronic 
microdevices. 
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35. The melhod of claim 29 ^vherein die arrayed data is 
generated from a device chosen from the group comprising: hybridized 
spotted cDNA microarrays, synthesized oligonucleotide arrays, spotted 
oligonucleotide arrays, peptide nucleotide assays, single nucleotide 

5 polymorphism (SNP) arrays, carbohydrate arrays, glycoprotein arrays,, 
protein arrays, proteomic arrays, tissue arrays, antibody arrays, antigen 
arrays, bioassays, sequencing microarrays, sequencing by hybridization 
(SBH) microarrays, siSNA duplexes, RNAi arrays glass-based arrays, 
nylon membrane arrays, thin fihn arrays, polymer-substrate arrays, 
10 capillary electrophoresis arrays, genospectral arrays, electronic arrays, 
bead arrays, quantum dot arrays, glycan arrays, spotted wells, spotted well 
plates. 

36. A system for analyzing a signal pattem to identify events of 
interest within the signal pattem comprising: 

15 an expressor function input unit; 

a preconditioner unit; 
an active interferometric coupler; and 
a resonant marker detector. 

37. A S3^tem for analyzing an arrayed signal pattem generated 
20 by an arrayed platform device to identify events of interest within the 

signal pattem comprising: 

an expressor ftmction input unit for inputting expressor functions; 

a preconditioner unit for preconditioning the arrayed signal pattem 
so as to convert the arrayed signal pattem to a spectral domain in which 
25 spectral harmonics parameterize events of interest to a pre-determined 
dynamical system; 

a coupler unit for convolving the preconditioned signal pattem and 
the expressor frmctions so as to interferometrically enhance portions of the 
preconditioned signal pattem associated wiA events of interest, if any. 
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present^vilfain the preconditioned signal pattern; and 

a resonant marker detector for identifying occurrence of events of 
interest within fhe enhanced signal pattern. 

38. The systsem of claim 37 wherein the signal pattern is imaged 
5 and pixilated. 

39. ITie system.of claim 37 wherein the e?q>ressor functions are 
designed so as to be capable of extracting one or more of spatial^ spatio- 
temporal or spectral invariants of events of interest associated with the 
array device. 

1 0 40. A method for anal3^zing an arrayed signal pattern generated 

by -an array device to identify events of interest within the signal pattern 
comprising the steps of: 

inputting expressor fimctions capable of extracting spectral 
invariants of events of interest associated with the array device; 
15 preconditioning fhe arrayed signal pattern so as to convert the 

armyed signal pattern to a spectral domain; 

convolving the preconditioned signal pattern and the expressor 
functions so as to interferometrically enhance portions of the 
preconditioned signal pattern associated with events of interest, if any, 
20 present within the preconditioned signal pattern; and 

identifying events of interest within the enhanced signal pattern. 

41. The method of claim 40 wherein the step of convolving the 
preconditioned sigoal pattem and the expressor fimctions includes the step 
of performing a reverberation convergence between the preconditioned' 

25 signal pattem and the expressor functions to achieve a resonance state. 

42. The method of claim 41 wherein the step of performing a 
. reverberation convergence is performed using an open loop control process 

that terminates upon detection of a predetermined condition. 

43. The method of claim 40 wherein the predetermined 



wo 2004/102456 PCT/US2003/021525 

-49- 



condition is a resonant marker corresponding to an event of interest. 

44. The method of claim 40 wherein the step of convolving the 
preconditioned signal pattern and the expressor functions inchides the step 
of a destructive interference between the preconditioned signal pattern and 

5 the expressor functions. 

45. The method of claim 40 flie step of convolving the 
preconditioned signal pattern and the expressor functions includes Ihe step 
of a constructive interference between the preconditioned signal pattern 
and the expressor functions. 

10 46. A system for analyzing a signal pattern to identify events of 

interest within the signal pattern comprising: 
an expressor function input unit; and 

an adaptive interferometric coupler operative to perform signal 
preconditioning, convolution coupling and resonant marker detection using 
15 the expressor functions. 

47. A system for analyzing an arrayed signal pattern generated 
by an arrayed platform device to identify events of interest within the 
signal pattern comprising: 

an expressor function input unit for inputting expressor functions 
20 C£qpable of extracting spectral invariants of events of interest associated 
with tiie anrayed signal pattern; and 

an adaptive interferometric coupler for preconditioning the armyed 
signal pattern so as to convert the arrayed signal pattern to a spectral 
domain to a spectral domain in which spectral harmonics parameterize 
25 events of interest to a pre-determined dynamical system while convolving 
the signal pattern and the expressor functions so as to interferometrically 
enhance of portions of the signal pattern and to identify events of interest, 
if any, within the signal pattem. 

48. A method for analyzing an arrayed signal pattem generated , 
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by an arrayed platform device to identify events of interest within the 
signal pattern comprising the st^s of: 

inputting expressor functions cq>able of extracting spectral 
invariants of events of interest associated with tiie array device; and 

preconditioning the arrayed signal pattem so as to convert the 
arrayed signal pattem to a spectral domain to a spectral domain to a 
spectral domain in which spectral harmonics parameterize events of 
interest to a pre-determined dynamical system while convolving the signal 
pattem and the e^ressor functions so as to hiterferometrically enhance 
portions of the signal pattem and to identify events of interest within the 
signal pattem. 

49. The method of claim 48 wherein the step of preconditioning 
the signal pattem includes the step of performing a closed loop 
reverberation convergence between the signal pattem, as it is being 
preconditioned, and the introduced e^qpressor functions to achieve a 
resonance state, terminated upon detection of a predetermined condition. 

50. The method of claim 48 wherein the predetermined 
condition is a resonant marker. 

51. A system for analyzing a signal pattem to identify events of 
interest within the signal pattem comprising: 

^a preconditioner unit; and 

an exprcssoT function adaptation unit operative to perform 
convolution coupling and resonant marker detection using canonical 
expressor functions. 

52. A system for analyzing an arrayed signal pattem generated 
by an arrayed platform device to identify events of interest within the 
signal pattem comprising: 

a preconditioner for preconditioning the arrayed signal pattem so as 
to convert the arrayed signal pattem to a spectral domain to a spectral 
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domain in which spectral harmonics parameterize events of interest to a 
pre-determined dynamical system; and 

an e^ressor function adaptation unit for generating preconditioned 
expressor functions based on canonical expressor ftinctions and based on 
5 the preconditioned signal pattern while extracting spectral mvariants of 
events of interest associated with the array device and interferometrically 
enhancing portions the preconditioned signal pattern so as to identify 
events of interest, if any, within the signal pattern. 

53. A method for analyzing an arrayed signal pattern generated 
10 by an array device to identify events of interest within the signal pattern 

comprising the steps of: 

preconditioning the arrayed signal pattern so as to convert the 
arrayed signal pattern to a spectral domain to a spectral domain in which 
spectral harmonics parameterize events of interest to a pre-determined 

1 5 dynamical system; and 

generating preconditioned expressor functions based on canonical 
e3q}ressor functions and based on the preconditioned signal pattern while 
extracting spectral invariants of events of -interest associated with the array 
device and interferometrically enhancing portions the preconditioned 

20 signal pattern to identify events of mterest, if any, wifliin the signal pattern. 

54. The method of claim 53 wherein the step of generating 
preconditioned e7q)ressor functions is performed so as to implicitly achieve 
a reverberation convergence between the preconditioned signal pattern and 
the expressor functions to achieve a resonance state. 

25 55. The method of claim 53 wherein the step of generating 

preconditioned expressor functions is performed using a closed loop 
process. 

56. A system for analyzing a signal pattern to identify events of 
interest within tiie signal pattem comprising: 
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an expressor function input unit; 

a preconditioner unit; 

an iterative interferometric coupler; and 

an adaptive controller operative to. control the iterative 
5 interferometric coupler to perform convolution coupling and resonant 
maricer detection using the input expressor fimctions, 

57. A system for analyzing an arrayed signal pattern generated 
by an arrayed platform device to identify events of interest within the 
signal pattern comprising: 

10 an expressor function input unit for inputting expressor functions 

capable of extracting spectral invariants of events of interest associated 
with the array device; 

a preconditioner for preconditioning the arrayed signal pattern so as 
to convert the arrayed signal pattern to a spectral domain to a spectral 

1 5 domain in which spectral harmonics parameterize events of interest to a 
pre-determined dynamical system; 

an iterative interferometric coupler for convolving the 
preconditioned signal pattern and the expressor functions so as to 
interferometrically enhance portions of the preconditioned signal pattern 

20 associated with events of interest, if any, present within the preconditioned 
signal pattern; and 

an adaptive controller for controlling the coupler to iteratively and 
selectively convolve expressor functions to the preconditioned signal 
pattern imtil a predetermined degree of convergence is achieved so as to 

25 identifying events of interest within the enhanced signal pattern. 

58. A method for analyzing an arrayed signal pattern generated 
by an arrayed platform device to identify events of interest within the 
signal pattern comprising the steps of: 

inputting expressor fimctions capable of extracting spectral 
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invariants of events of interest associated with the array device; 

preconditioning the arrayed signal pattern so as to convert the 
arrayed signed pattern to a spectral domain; and 

iteratively and selectively convolving the preconditioned signal 
5 pattern and the expressor functions so as to interferometrically enhance 
portions of the preconditioned signal pattern associated with events of 
interest, if any, present within the preconditioned signal pattern until a 
predetermined degree of convergence is achieved so as to identifying 
events of interest within the enhanced signal pattern. 

1 0 59. The method of claim 58 wherein the step of iteratively and 

selectively convolving the preconditioned signal pattern and the expressor 
functions is performed by repeatedly convolving a single fixed e^ressor 
function to the preconditioned signal pattern using reverberant 
convergence imtil the predetermined degree of convergence is achieved. 

15 60. The method of claim 58 wherein the step of iteratively and 

selectively convolving the preconditioned signal pattern and the expressor 
functions is performed by selectively convolving a set of different fixed 
expressor fimctions to the preconditioned signal pattern using reverberant 
convergence until the predetermined degree of convergence is achieved. 

20 61. The method of claim 58 wherein the step of iteratively and 

selectively convolving the preconditioned si^al pattern and the expressor 
functions is performed by selectively modifying and then convolving a 
dynamic expressor function to the preconditioned signal pattern using 
reverberant convergence until the predetermined degree of convergence is 

25 achieved. 

62. A system for analyzing an arrayed signal pattern generated 
by an arrayed platform device to identify events of interest within the 
signal pattern comprising: 

means for inputting e^ressor fimctions capable of extracting 
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spectral invariants of events of interest associated witihi the array device; 

means for preconditioning the arrayed signal pattern so as to 
convert the arrayed signal pattern to a spectral domain; and 

means for iteratively and selectively convolving the preconditioned 
5 signal pattern and the expressor functions so as to interferometrically 
enhance portions of the preconditioned signal pattern associated with 
events of interest, if any, present within the preconditioned signal pattern 
imtil a predetermined degree of convergence is achieved so as to 
identifying events of interest within the enhanced signal pattern. 
10 63. A computer code product for actively analyzing a signal 

pattern representative of arrayed data to identify events of interest therein, 
comprising: 

computer code that inputs a signal pattern representative of arrayed 

data; 

15 computer code that generates resonance pattems based on 

interference between syntiiietic noise and the signal pattern; and 

computer code that detects resonances within the resonance 
pattems associated with events of interest 

64. A computer code product for analyzing an arrayed signal 
20 pattern generated by an arrayed platform device to identify events of 
interest within the signal pattern, comprising: 

computer code that preconditions the arrayed signal pattem based 
on input expressor &nctions so as to convert the arrayed signal pattem to a 
spectral domain in which spectral harmonics parameterize events of 
25 interest to a pre-determined dynamical system; 

computer code that convolves the preconditioned signal pattem and 
the e3q>ressor fimctions so as to interferometrically enhance portions of the 
preconditioned signal pattem associated with events of interest, if any, 
present within the preconditioned signal pattem; and 
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computer code that identifies the occurrence of events of interest 
within the enhanced signal pattern. 
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